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Abstract 

There is a proposition that the possible occurrence of electrostatic discharges inside an oxygen regulator 
between charged insulating components, such as diaphragms, and earthed metal components, such as the valve 
stem, could be a potential candidate responsible for oxygen regulator fires. In this paper, the electric field which 
might be produced at the tip of the valve stem due to charge produced on the diaphragm of a typical oxygen 
regulator by deposition of charged dust particles or by their impact on the diaphragm or both has been evaluated. 
The diaphragm, in many cases, is backed partially by a concentric earthed metal disk. The diaphragm has been 
modelled as an insulating disk. The tip of the valve stem has been modelled as an earthed sphere which is very 
small compared to the insulating disk and the backing earthed metal disk. An analytical expression for the 
electric field produced at the point of the sphere nearest to the charged surface of the insulator has been derived. 
Our expression takes into account not only the effect of the charge on the insulating disk, but also that due to the 
presence of the backing earthed metal disk. Results for the magnitude of this electric field have been computed 
for the case of a sample oxygen regulator. An expression for the critical charge density on the surface of the 
insulating disk has been derived, and its value has been obtained for the case of our sample oxygen regulator. 
The electric field inside the non-metal diaphragm has also been evaluated. It appears from our analyses that 
there is a possibility that an electrostatic discharge might occur inside an oxygen regulator, and with an 
enriched-oxygen atmosphere being present there, such a discharge could also lead to a fire incident. 


1. Introduction 

Oxygen regulators, used in places that require oxygen gas of high purity for their operations, 
such as, welding shops, steel plants and medical facilities, have been responsible for many 
fire accidents. A number of such incidents have been reported [1 - 5], In some of these 
incidents again, the mechanism of ignition could not be firmly established [4, 5], This 
indicates that apart from the well-known mechanisms mostly used to explain oxygen 
regulator fires, namely, ‘particle impact’ [1, 6], ‘mechanical impact among regulator 
components’ [6] and ‘adiabatic compression’ [6], there could be other mechanisms which 
might also be responsible for such incidents. Two possible mechanisms have been suggested. 
These are ‘flow friction’ [6 - 8] and ‘tribocharging of particle contaminants’ [9, 10]. 
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The last mechanism, i.e., ‘tribocharging of particle contaminants’ [9, 10], suggests 
that particle contaminants, while flowing through the oxygen system, might get 
electrostatically charged due to contact with the inner walls of the system and due to impact 
amongst themselves. These charged particles could eventually accumulate on the diaphragm 
of the regulator. This might result in a brush discharge when, during opening of the regulator, 
the diaphragm comes very close to the valve. The discharge might eventually cause fire due 
to ignition of dispersed dust particles and of non-metallic parts of the regulator present in the 
enriched-oxygen atmosphere of the regulator. 

An electrostatic discharge inside a regulator is also possible if the insulating 
diaphragm itself gets charged by the impact of dust particles. That insulating polymer 
surfaces get charged on being impacted by particles has been experimentally shown using 
five different polymers (fibreglass, Lucite, Rulon J, Lexan and Teflon) and Mars simulant 
particles having diameters less than 5 pm by making the particles travel at 20 m/s in an 
atmosphere of CCf at 10 mbar pressure [11]. 


Whatever the mechanism of charge generation on the diaphragm surface, whether the 
settling of charged particles on its surface or the charging of the surface itself by impact of 
particles or both, we have a situation where there is a charge distribution on the surface of an 
insulator backed partially by an earthed metal plate which charged-insulator-backed-by- 
earthed-metal-plate system could accidentally come close enough to an earthed electrode 
(which in this case is the valve stem) to cause an electrostatic discharge. Hence, it is 
important to know the electric field that could be produced in such situations. Since the 
magnitude of the electric field is maximum at the tip of the electrode, we shall be interested 
in this magnitude as a discharge will ensue if it exceeds the breakdown strength of the 
medium. Heidelberg [12] has shown that the magnitude of the electric field at the lowest 
point of an earthed conducting sphere placed above a uniformly charged circular insulating 
surface oriented such that the straight line joining its centre with that of the sphere is 
perpendicular to the charged surface at its centre is given by 
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where a is the density of charge on the insulating surface, R is the radius of the insulating 
surface, a is the radius of the earthed sphere, D is the distance between the centre of the 
sphere and that of the charged surface and s is the permittivity of the medium. He has also 
stated that if the insulating disk is backed by an earthed metal plate of the same radius and if 
the distance ;/ between the charged insulating surface and the metal surface is small compared 
to the distance between sphere and the insulator, then mathematically, the case may be treated 
approximately by considering the difference between the effects of two single layers of 
charge of same magnitude on the sphere, so that the electric field at the lowest point of the 
earthed sphere is given by 
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Ev = Y e [/(a, D, R) - f(a, D + rj.R )] ... (B), 
where the function/has been defined in equation (A). 

Oza et al [13] have shown by Computational Fluid Dynamics (CFD) simulation that 
the concentration of dust particles inside the regulator becomes more or less uniform after its 
closure. These dust particles gradually settle on the diaphragm surface. Thus, the distribution 
of charge on the diaphragm surface is more or less uniform. 

However, in many cases, the diaphragm is only partially backed by a metal disk as 
shown in Fig. 1 so that spring-loading is possible. 




Fig. 1. Photograph of a Regulator Diaphragm Showing the Backing Metal Plate 

There is also a metal disk in front of the diaphragm by which the valve stem can be pushed. 
For ease of analysis, we have ignored its presence as its purpose being only to push the valve 
stem, its radius can be brought down to that of the valve stem which is very small compared 
to both the diaphragm and the backing earthed metal disk. Due to the presence of this backing 
plate - which often gets earthed when a part of an oxygen system - covering some area at the 
centre of the diaphragm and because of the dust particles being charged, it may be expected 
that dust particles will settle in greater proportions at and around the centre of the diaphragm 
than towards its periphery. This is because the electric lines of force resulting from the 
charged particles will tend to get more concentrated towards the earthed backing plate, 
thereby causing the particles to be drawn towards it. Hence, the density of charge on the 
diaphragm surface may be expected to be the maximum at the centre and gradually reducing 
with increasing distance from it. In this work, we investigate the electric field produced at the 
tip of the valve stem due to this nature of charge distribution on the diaphragm surface. 

Till now, mathematical analysis of the electric field produced due to a charged 
insulating plate has been done only for cases pertaining to a uniform charge distribution. 
These results cannot be directly used in our problem as the charge distribution on the 
insulator (the diaphragm) for our case is non-uniform. We, in this work, mathematically 
analyse the electric field for a charge distribution which is such that the density of charge is 
the maximum at the centre of the insulator and gradually diminishing with increasing 
distance from the centre. Such a charge distribution may be faithfully represented 
mathematically by the Gaussian function. We model the tip of the valve stem as a very small 
earthed sphere, the diaphragm as an insulating disk and the earthed backing plate as an 
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earthed metal disk smaller than the insulating disk and placed concentrically with the 
insulating disk. 

2. Statement of the Mathematical Problem 
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Fig. 2. An Earthed Sphere near a Charged Insulator Partially Backed by an Earthed Metal Disk 

We consider a system consisting of an earthed conducting sphere placed near a positively 
charged insulating disk of thickness d, which is partially backed by a concentric earthed metal 
disk, such that the straight line joining the centres of the sphere and the insulator is 
perpendicular to the insulating surface at its centre as shown in Fig. 2. Let a represent the 
radius of the sphere, R the radius of the insulating surface, b {b <R) the radius of the backing 
earthed metal disk, D the distance between the centres of the sphere and the insulator, and 
er(r, 9) the surface density of charge on the insulating surface. It is required to find the 
magnitude of the electric field Ep at the lowest point P of the sphere for the surface charge 
density on the insulator given by 



where cr 0 is the maximum density of charge on the insulator and s is a parameter which 
controls the width of the “bell” of this bell-shaped charge distribution. 


3. Solution 


First, let us find out the electric field at P due to the charge on the insulator without 
considering the effect of the backing earthed metal plate. We shall use the method of images 
[14, 15] for our analysis. The Cartesian coordinate axes are chosen as shown in Fig. 2. Let us 
consider an element of charge dq at a point X as shown in Fig. 2. Let OiX = r, and the angle 
subtended by OiX with the ;c-axis be 9. Then, we have as indicated in Fig. 2 
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dq — <j(r, 9) ■ rdO ■ dr — a(r, 9)rdrdO 

Let the image charge due to the element of charge dq be located at I. Then, we have [14, 15] 
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Let Z be any point on the z-axis. Then the potential at Z due to the element of charge dq and 
its image dq' is given by 
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Again, we have 
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From equations (4) and (5), we have 
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The component of the electric field at Z along the z-axis is given by 
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Now, we have 
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Putting z = (D-a ) in equation (10), we have the electric field at the point P (the lowest point 
of the earthed sphere) due to the element of charge dq and its image dq' as 
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Hence, the net electric field at point P due the entire charge on the insulating surface is given 
by 
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Using the expression for er(r, 0) given by equation (1) in equation (12), we have 
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Equations (14) and (15) give expressions for the electric field at P, the lowest point of the 
earthed sphere (as shown in Fig. 2), due only to the charge on the surface of the insulator. 
Now, we shall proceed to find the net electric field at P taking into account the effect of the 
earthed metallic backing disk. No attempt has yet been made to find an exact analytical 
expression for this net field. We shall find here analytical expressions for the net electric field 
at P. 
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Finding the Net Electric Field at P Taking into Account the Effect of the Earthed Metallic 
Backing Disk 

As has already been mentioned in section 1, Heidelberg [12] has stated that for the case of a 
charged insulating disk being backed by an earthed metal disk of the same radius, if the 
distance between the charged insulating surface and the backing metal surface is small 
compared to the that between sphere and the insulator, then the case may be treated 
approximately by considering the difference between the effects of two single layers of 
charge of same magnitude on the sphere, so that the electric field at the lowest point of the 
earthed sphere is given by equation (B) with the function / appearing in that equation being 
defined by equation (A). This may be done because the charges on the charged surface of the 
insulating disk continue to induce charges in the earthed metal disk until the electric field 
inside the earthed metal disk is zero and the electric field on the surface of the earthed metal 
disk is perpendicular to this surface, and this state of the electric field is achieved only if a 
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charge of a surface density which is same in magnitude as, but opposite in sign to that of the 
charge on the charged surface of the insulating disk is induced on the surface of the earthed 
metal disk facing the charged surface of the insulating disk and the charge induced on the 
other bounding surfaces of the earthed metal disk is zero. (This is so only if the distance 
between the surface of charge on the charged surface of the insulating disk and the surface of 
the earthed metal disk facing this surface of charge is very small compared to the radii of 
these two surfaces.) The charge induced on the surface of the backing earthed metal facing 
the charged surface of the insulating disk is too tightly bound to the charge on the insulating 
disk to be able to move freely in the metal, and hence, it may be treated as though it were 
residing on an insulating surface. Using this approach and considering the charge induced on 
the earthed backing disk to be almost solely due to the charge residing on the insulator 
immediately above the metal disk (and considering the effect of the charge on the insulator 
residing beyond the distance b from the centre to be negligible so far as the production of the 
induced charge is concerned), an assumption which might be considered valid to a good 
approximation for S « D. R. b, which is true for our case, the net electric field at P is given 
by 

E Pnet =-^-[<t>(D,a,R,s)-</>(D + S,a,b,s)\ ... (16) 

2m 

or,E Pnet =-^[(//(g,a,R,s)-(//(g + S,a,b,s)\ ... (17) 

2sa 

If we wish to consider the image charges produced in the earthed metal disk instead of the 
charges induced on its upper surface - the surface facing the charges on the insulating disk - 
to calculate the electric field, we will have the following expression for the net electric field 
at P: 

E p n et = -^-[</>(D,a,R,s)-t/>(D + 2S,a,b,s)\ ... (18) 

2m 

or,E Pnet = -^[i//(g,a,R,s)-i//(g+2S,a,b,s)\ ... (19) 

2m 

This is because each element of charge (having a charge o(r, 0).rdrd0) on the charged surface 
of the insulating disk upto the distance b from the centre produces an image charge which has 
a charge of a magnitude same as, and a sign opposite to that of the element of charge itself on 
the side of the earthed metal disk surface facing the charged surface of the insulating disk 
opposite to that of the charged surface of the insulating disk at a distance S from this earthed 
metal disk surface on the straight line which passes through the said element of charge and is 
perpendicular to this surface of the earthed metal disk [14, 15]. We note that for the elements 
of charge on the insulating disk (each having a charge a(r, 9).rdrd0) that lie towards the 
periphery of the earthed metal disk, the surface of the earthed metal disk facing these charges 
cannot be considered to be really forming a closed surface which bound spaces internally and 
externally, which is a condition which is necessary to ensure that the electric field produced 
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on this metal surface due to each of these elements of charge is the same as that obtained on 
this surface by adding the electric fields due to this element of charge and another charge 
which we call its image which has a magnitude and a sign of charge and a position such that 
the surface described by this surface of the earthed metal disk (the surface of the earthed 
metal disk facing the charged surface of the insulating disk) is an equipotential surface and 
has a potential same as that of the metal disk (which in this case is zero). This requirement in 
using the method of images is demanded by the following theorem: 

“If the value of the potential V is known at every point on a number of closed surfaces by 
which a space is bounded internally and externally, there is only one value for V at every 
point of this intervening space, which satisfies the condition that V l V either vanishes or has 
an assigned value at every point of this space.” [15] 

Thus, in using the method of images for the elements of charge on the charged surface of the 
insulating disk that lie towards the periphery of the earthed metal disk, we have made the 
assumption that the surface of the earthed metal disk facing these elements of charge to be 
forming a closed surface which bound spaces internally and externally despite these elements 
of charge being near to the boundary of this surface. This assumption is valid as each of these 
elements of charge has infinitesimal charge and, hence, has a negligible effect on the points 
in the space on the side of the earthed metal disk opposite to that in which each of these 
elements of charge is located. These images, again, give rise to the distribution profile of the 
induced charge in the earthed metal disk as has been described above— each point on the 
surface of the earthed metal disk facing the charged surface of the insulating disk having a 
surface charge density having magnitude same as and sign opposite to that at the point on the 
charged surface of the insulating disk right above it and the other surfaces of the earthed 
metal disk bounding it having zero charge. 

We assume the effect of the charge on the insulator beyond the distance b from the 
centre in producing image charges in the earthed backing disk to be negligible. The 
assumption is valid to a good approximation in our case since we have here S « R , b. That it 
is true may be understood by reasoning that in such a case, the charges on the surface of the 
insulator right above the earthed metal disk will have the maximum effect in producing the 
induced charges on the surface of the earthed metal disk and the charges on the surface of the 
insulator beyond the distance b from the centre will have little effect in producing these 
induced charges, and that image charges have the effect of only producing this induced 
charge, and hence, it must be so that image charges are produced in the earthed metal disk 
almost solely due to the charges on the surface of the insulator right above the earthed metal 
disk, and very little due to the charges on the surface of the insulator beyond the distance b 
from the centre. Our approximation— that of neglecting the effect of the charges on the 
surface of the insulator beyond distance b from the centre in producing the image charges in 
the earthed metal disk— can only have the effect of varying, if at all, only slightly the electric 
field and the induced-charge distribution profile on the surface of the earthed metal disk, but 
these variations can cause very little changes to the electric field at P— the question that we 
are interested in here. 
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In the above analysis, the effect of reflection of the image charges produced in the 
earthed sphere on the earthed backing disk and the reflection of these new images back on the 
sphere, and such reflection of image charges going on on the two earthed metal surfaces ad 
infinitum has been neglected. Now, we shall try to take care of this effect. We do our analysis 
below. 

Our system may be reproduced as follows: 



Fig. 3. An Earthed Sphere in Front of a Charged Insulating Disk Partially Backed by an Earthed Metal Disk 
(Figure Depicting the Longitudinal Section of the System Shown in Fig. 2) 

As has been shown in Fig. 3, we assume that a, S « D , b, R. The assumption that a, S « D 
is valid for our system of interest for all cases except when the sphere is very close to the 
insulating disk. The other assumptions are always valid for our system. We draw the 
following figure (Fig. 4) for ease of visualization in doing our analysis. 
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Fig. 4. Figure Meant for Ease of Visualization in Using the Method of Images such that the Effect of Reflection 
of the Image Charges Produced in the Earthed Metal Sphere on the Earthed Metal Backing Disk and the 
Reflection of These New Image Charges Back on the Sphere, and Such Reflection of Image Charges Going On 
on The Two Earthed Metal Surfaces Ad Infinitum Is Taken Care Of (A Denoting the Earthed Metal Sphere in 
the System of Our Interest, B Denoting the Surface of Charge on the Insulating Disk, B' Denoting the Surface of 
Image Charges Produced Due To Reflection of The Charges on the Surface of Charge B Upto the Distance b 
From the Centre on the Earthed Metal Disk, /, and i 2 Denoting the Angles Subtended by the Surfaces B and B' 
Respectively at Centre O of the Earthed Metal Sphere A) 

As has been explained above, B' will have the same charge distribution profile (in 
magnitude) as B has up to the distance b from the centre, but the charges on B' will differ 
from those on B in sign. 

The electric field at point P will be produced by the charge surfaces B and B', the 
image charges q f and <y, produced by these in sphere A, the images produced in the earthed 

metal disk by these image charges, the images produced by these back in sphere A, their 
images produced again in the earthed metal disk, and so going on ad infinitum. The electric 
field at P is, thus, given by 


^P,net ^B,q h 


+ E„' + E 


A-la 


+ E f , ...(20), 


where 


E B q is the electric field produced at P by the charge surface B and its image charge in 
sphere A, 

E B , q is the electric field produced at P by the charge surface B' and its image charge q h in 
sphere A, 

E, j e i is the electric field produced at P by all image charges in sphere A except 
G,(= <?,, +q h ), and 
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Ej is the electric field produced at P by all image charges at point O'. 
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We have, thus, found expressions for the first two terms for E Pnet given in equation (20). 
Now, we shall proceed to find out expressions for the remaining two terms given in this 
equation, namely, E,^ _/ g| j and E t . In order that we may do so, we start by finding out 

expressions for q t and q h . 

We first proceed to find out q h , i.e., the image charge produced by the charge surface 
B in sphere A. Using equation (3), we have 


= - 


a 


a/ D 2 + r 2 


(j(r,0)rdrdG ... (23) 


This image charge is located on the line segment joining the centre O of sphere A and the 


element of charge dq on the surface B at a distance 


a~ 


a/z) 


2 . 2 

+ r 


from O. Since we have 


assumed that a « D , R , b, and also that a itself is very small, this image charge may be 
assumed to be located almost at the centre O of the sphere. The total image charge produced 
in A by the charge surface B is given by 


2 7tR 
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a{r,6)rdrd0 ... (24) 
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Similarly, it can be shown that the image charge produced in sphere A by the surface of 
charge B' is given by 


(D+2S) 2 

q t = y[2ascr 0 7ve 2,2 



1 (Z) + 2£) 2 " 
2’ 2? y 


1 b 2 + {D + 2Sf\\ 

2’ 2? Jj 


...(27) 


and that this charge too may be assumed to be located almost at the centre O of this sphere. 

We, thus, have a charge Q t (= c/. + q t ) located almost at the centre O of the earthed 

metal sphere A. O is located at a distance (D + 6) from the earthed metal disk. The charge Q\ 
will reflect on the surface of the earthed metal disk, and will, thus, produce an image charge 
of value (-<2i) at point O' which is located at a distance (D + 6) from the earthed metal disk 
on the side of it opposite to that where sphere A, and hence, point O and charge Q i are 
located. (Kindly refer to fig. 4.) This image charge of value (- Q \) produced at O' will reflect 
on the surface of the earthed sphere A, and will, thus, produce in it an image charge of value 

a a 2 

- O , at a point located at a distance - from O on the line joining O and O'. 

2 (D + S) 2 (D + S) 

Since we have assumed for our system a « D , R , b, and also that a itself is very small, the 
distance of this image charge from O may be considered negligible, and this image charge 

may be assumed to be located almost at O. This image charge of value ^ ^ — Q l located at 

O will reflect on the earthed metal disk, and will, thus, produce at O' an image charge of 

value j- — Qi | • Reflections of image charges will continue to occur on the earthed 

metal sphere A and on the earthed metal disk following the manner described here ad 
infinitum. Since we have assumed for our system that a « D , R , b , and also that a itself is 
very small, we consider for all image charges in the earthed sphere A their distances from the 
sphere’s centre O to be negligible. In other words, we assume in our calculations that in the 
earthed sphere A, all the images are produced almost at its centre O. This also results in 
underestimation of the magnitudes of the image charges by a certain extent. However, 
because of our assumption that a « D , R , b, we may consider this underestimation to be 
negligible. Based on this assumption, we have an infinite number of image charges produced 
at points O and O' as shown below in Table I. 
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Table I: Image Charges Produced at Points O and O' 


Image Charges Produced at Point 0 

Image Charges Produced at Point 0 ' 
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[image produced due to reflection of Q x on the 
earthed metal disk] 
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[image produced due to reflection of (- Q \) on the 
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on the earthed metal disk] 

• 

• 


We then have 


E Ia~{Qi 


4 7iar 


• 0 ,+ 


0i+ ■ 


2(D + S) 1 {2 {D + S)f ‘ {2(D + S)\ 


13 ^1 


0 , 


{2 (D + S)) 


4 VSl 


01 H-tO oo terms 


... (28) 


°r,E, 


fa 


0 , 


4/rarr 




2(D + f>) + {2(£> +f>)} 2 ^ {2(D + <?)} 3 ^ {2(D + <?)}' 


- + 


a 


+ ■ 


- H— to oo terms 


19 




Gi 


4/rar 


2 (P + S) 


1- 


a 


2{D + S) 


a 


0, 


or ’ £ /,-{a} = - 


or ’ E u-{Q4 


4;ra7 2 |2(.D + c>)-aj 

0. 

2(D+S) 

a 

4^a7 2 {2(D + S)-ci} 

1 

a 

1 4xai 2 

{2{D + S)-a} 

1 

a 


Gi 


... (29) 


4^a7 2 {2(Z) + c>)-a} 


- 42as(7 ri 7ie 2s2 \ r 


Z) 2A 


v 2 25' y 


-r 


r i r 2 +d 2A 


v 2 2r , 


■ a/2, 


(£>+2<S) z 


+ v las cr 0 7ie 


1 (D + 25) 
2’ 2? 


,2 A A 

-r 


1 6 2 + (£> + 2£) 

2 ?" 


,2 A 


V ““ J V 

0! = q h + q t , and substituting expressions for q t and <y, given by equations (26) and (27) 
respectively] 

r 1 a 

or ' E i A -{Q4 = 


4nm 2 {2(D + S)- 


42a serene 


a] 

(D+2S) Z 


2s* 


1 (D + 2Sf) ( 1 b 2 + {D + 28) 


v 2 2s l j 


-r 


a 4 


2s 2 


- V2 as<j n Ke lsl \ T 


f \_ 

2 ’ 2 ? 


-r 


1 R- + D 2 

v 2 25 y 


or ’^-{a} = 


1 Cl nr 

'-)’(/ \ v * v 2^asGrJi 

4 nsa 1 {2(D+ <?)-. 1 


(D+2<?r 


1 (D + 23) 


v 2 2s" y 


,2A f 

-r 


1 b 2 +(D + 2S) 
2 ’ 


2s* 


-e 2 ’ 2 r 


f \_ 

2 ’ 2 ? 


-r 


1 R 2 + D 2 

2’ 2 s 2 


or, E, 


-e^ 2 r 


_ &0 1 

f 5 ^ a 

(D+25 ) 2 

2j2 Jr 

f 1 (D + 2S) 2 'l 

_r 

f 1 6 2 + (D + 2^) 2 Y| 

lal 2j2e 1 

\a) {2(p + d)-a) 

r r 

^2’ 2s 2 j 


l 2 ’ 2 * 2 Jj 


I ^ 

v 2’25 y 


-r 


l R +D 2 

v 2’ 25 y 


20 



or, £ 




a o 


a 

(D+2S) 2 

^ 2s2 i r 

'i (d + 2 s ) 2 N 

_r 

7 (£> + 27) 2 +Zr7 

2V2f 

UJ 

{2 (D + S)-a} 

1 

2’ 2s 2 

1 

l 2 ’ 2 * 2 Jj 


-e 2 ’ 2 T 


... (30) 


E, =■ 


O Z)_ 2A 

v 2 ’2? j 


-r 


f i d 2 + r 2A 


v 2 25 2 y 


4^{2(Z) + S) - a } 9 
-to oo terms] 

r 61 

or, E, =- 


-0 _«_n_ 


a- 


2(Z) + O’) 1 {2(£> + £)} 2 1 (2(D + ^)} 


■fii 


...(31) 


4^-£-{2(d + ^)-o} 2 
0, 


, a a a 

1 + 7 -r + —;- ttt + -r-p-— to oo terms 


2 (D + S) {2 (D + S)} 2 {2{D + S)f 

1 


or,E, =- 


or, E, =- 


4^o{2(D + 7) - a }° | 

a 

1 2 (D + (7)} 

2 (D + S) 

4;ro{2(.D + c>)-o}' {2(D + 7 )- o} 

1 

2(D + S) 

4^{2(D + <7)-o} 2 

_{2(d + 8)-ci} 

1 

2(D + S) 


Qi 


... (32) 


4tts{2{d + S)-of { 2{D + S)-a) 


- -Jlascjrjze'-’ 1 \ T 


( 1 r\- \ f | _|_ 2) 2 ^ 


I D 

2'2s 2 


-H 


2 2s 2 


■JL 


( p+isr 


+ V las o^Tie 


{r 

7 (D + 2S) 2 ' 

-T 

1 

2’ 2 s 2 

V 


1 b 2 + (D + 2Sj 

2s 2 




[v 0! = q, t + q, , and substituting expressions for q t and q t given by equations (26) and (27) 
respectively] 


or, E, = - 


1 


2(D + 8 ) 


4tts{2(D + <?)- of ' {2 (D + S)- r 1 ' ^ asa * n 


(D+2Sr 


1 (D + 2d') 


v 2 2.S y 


a 

'V 


1 Zr+(£> + 2£) : 

2’ 27" 


-e 2 - 2 r 


I D 

2 ’ 2 ? 


-r 


1_ £ 2 + D 2 

y2’ 2s y 


21 



or, E, = — 


^0 I 

fs\ 

a 

2 

2 (D + S) 

242s 1 

yaj 

{2 (D + S)-a}_ 


_{2 {D + 5)-a} 


1 (D + 28) 


2s 2 


f 

-r 


y 


1 (£> + 2<V) 2 +lr 
2 ’ 2 ? 


... ( 33 ) 


-e 2 ’ 2 AT 


f \_ 

v 2’2V j 


-r 


A i d 2 + r 2 ^ 


v 2 2* , 


We now have expressions for all the terms of E p net given by equation (20) which equation is 
reproduced below for convenience. 

Ep.net = E riq . + E b , ^ + E Ia _ {a} + E, . . . (20), 


where E Bq is given by equation (21), E B , q is given by equation (22), E Ia _^ is given by 
equation (29) and by equation (30), and E { is given by equation (32) and by equation (33). 
These equations are reproduced below for convenience. 
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We have, thus, found our expression (by far the most accurate) for the electric field produced 
at the point P, the point of the earthed sphere in the system of our interest (the system 
constituted by an insulating disk of which one surface has charges on it and the surface 
opposite to it has on it an earthed metal disk and an earthed metal sphere approaching the 
surface of the disk on which charges are present as shown in fig. 2) which point is nearest to 
the charge surface on the insulating disk for a charge distribution profile on the surface of the 
insulating disk which has charges on it given by 

<y{r,0) = cr Q e ^ ... (1), 

where cr 0 is the maximum density of charge on the insulator and s is a parameter which 

controls the width of the “bell” of this bell-shaped charge distribution. We write this 
expression for convenience as follows. 
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4. To Find the Minimum Value of cr 0 Required for a Discharge to Occur at the Point of 
the Earthed Metal Sphere Nearest to the Charged Surface of the Insulating Disk 

To find the minimum value of cr 0 required for a discharge to occur at the point of the earthed 
metal sphere nearest to the charged surface of the insulating disk, we need to put the value of 
the breakdown field for E Pnet and also put D = a in equation (35). Then, we have 
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or, cr 0 = 2sE hreakdown [f 2 (8, a, A 1 . /?,.sj] '... (42) 


Equation (42) used with equation (41), thus, gives us the minimum value of cr 0 required for a 
discharge to occur. 

We may also find the minimum value of <j 0 required for a discharge to occur at the 

point of the earthed metal sphere nearest to the charged surface of the insulating disk when its 
centre is at a distance D from the charged surface of the insulating disk. This value can be 
found by putting the value of the breakdown electric field in equation (36). On doing so, we 
have 


tfolmin (D,8, a,R,b,s ) 


2sE breakdown 

A (D,8, a,R,b,s ) 


(42 A), 


where A iP, 8, a, R, b, s) is given by equation (37). 
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5. To Find the Electric Field at a Point Inside the Insulating Disk on the z-axis: 

We observe that the electric field inside the insulating disk is maximum at points on the z- 
axis as these points are nearest to the centre of the charged surface of the insulating disk 
where the density of charge is maximum as this centre is the origin of our Cartesian 
coordinate system. Hence, we find here the electric field at a point inside the insulating disk 
on the z-axis. We first find out the expression for the electric field at a point (0, 0, z) due only 
to the charge on the charged surface of the insulator, i.e., to the surface of charge B shown in 
Fig. 4. The x, y and z components of the electric field at a point (0, 0, z.) due to an element of 
charge at a point (r, 9) having a charge cr(r, 6).rdrd6 on the surface of charge B are as 
follows: 
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(45). 


Therefore, the v, y and z components of the electric field at the point (0, 0, z) due to the 
surface of charge B are as follows: 
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We have for two functions « and v 
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d(uv ) = udv + vdu 


:.uv=judv + jvdu 
,\udv = uv-\vdu 

_ _ 3 

Let u = e w and dv — w 2 dw 
Then, we have 
du — —e~ w dw,and 
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Then, we have 
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From equations (47), (49) and (51), it is evident that the electric field at a point (0, 0, z) due to 
the surface of charge B is given by 


o Q (z z 

Eb ~ 2c V|z| VR 2 + z 2 * “ 




2' 2s 2 


and it is directed along the z-axis. 
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It may be found by a way of analysis similar as the way shown above that the electric field at 
a point (0, 0, z) due to the surface of charge B' is given by 


E R ' = -■ 


On 


(z + 25) 


(z + 25) 


2£|\|(z + 25)| j h 2 + ( z + 25) 2 


■ e 2 s 2 


(z + 25) (z +2 s) 2 

Vis 


■ e 2 s 2 


r 


'1 (z + 25) 2 
.2' 2s 1 


- r 


'1 b 2 + (z + 25) 2 
2’ Is - 2 


(53) 


We have for any point inside the insulating disk on the /-axis, z = — |z| and 0 < |z| < 5. 
Hence, the electric fields at such a point due to the surfaces of charge B and B' are 
respectively given by 
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(54), and 
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Equations (54) and (55) may be expressed as shown below. 
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(57) 


We now proceed to find the electric field at a point (0, 0, z.) due to the image charges 
produced in the earthed metal sphere A, the images produced in the earthed metal disk by 
these image charges, the images produced by these image charges back in sphere A, their 
images produced again in the earthed metal disk, and so on going on ad infinitum. From table 
1, we have the following. 
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Total charge of all the image charges produced at 0 
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Qi ••• (59) [by evaluating as shown above] 


Ci = Qi, + qi 2 

Substituting here the expressions for q ti and q iz from equations (26) and (27) respectively, 
we have 
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Putting this expression for Q\ in equations (58) and (59), we get 
total charge of all the image charges produced at 0 
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The electric field at point (0, 0, z) due to all the image charges produced at O is given by 
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(63) [using equation (60)] 


The electric field at point (0, 0, z) due to all the image charges produced at O' is given by 
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(65) [using equation (61)] 


Therefore, the electric fields at any point inside the insulating disk, for which z = — |z| and 
0 < |z| < 5, due to all the image charges produced at O and due to all the image charges 
produced at O' are, respectively, given by 
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(67) [from equation (63)] 
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(69) [from equation (65)] 
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The net electric field at any point inside the insulating disk on the z-axis— for which 
z = — |z| and 0 < |z| < 8 — is given by 

^insulator,net ~ ^insulator,B ~b Einsulator, b' ~b Emulator,0 "b E insulator,O' (^0) 

where 
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6. Results and Discussions 

We have evaluated the electric field that could be produced at the tip of the valve stem of an 
actual oxygen regulator using our mathematical model. We present the results here. From 
measurements made for our sample oxygen regulator, we have obtained the following 
dimensions for our model: 

a = 0.8 mm, R = 27.75 mm, b = 17.4 mm and S = 1.8 mm. 

We have obtained results for three values of the maximum charge density cr 0 on the 

diaphragm surface: 0.2 x 10 11 C/mm 2 , 0.5 x 10 11 C/mm 2 and 1.5 x 10 11 C/mm 2 . We have 

tried with these values of charge density following the work of Heidelberg [16]. We have 

taken s = b = 17.4 mm in our calculations. In our calculations, we have used the electric 

permittivity of oxygen gas at NTP. The relative permittivity of oxygen at NTP is 1.00049 [6]. 

Multiplying it with the permittivity of free space (8.854 x 10" F/m), we obtain the electric 

12 

permittivity of oxygen gas at NTP as 8.858 x 10 F/m. 

First, we present (in Fig. 5) the results for the electric field produced at P without 
taking into account the effect of the earthed metallic backing disk, which we have obtained 
using equation (14). The graphs have been obtained using MATLAB R2014a software. 
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Fig. 5. Electric Field at P Obtained Without Taking into Account the Effect of the Earthed Metal Backing Plate 

In Fig. 5 as well as in Fig. 6, the values for the charge density given in the legends are the 
values of cr 0 . As expected, we find that the magnitude of the electric field at P is maximum 

when the sphere is almost in contact with the insulating disk, i.e., at D = 0.8 mm, and it 
decreases as the sphere goes farther and farther away from the disk, i.e., with increasing value 
of D. Also the magnitude of the electric field increases as the charge density on the insulating 
disk increases. 

Now, we present the results obtained for the net electric field at P discussed in section 

3. 


Net Electric Field With Backing Plate Approach-1 



Fig. 6. Results for Net Electric Field at P Obtained Using Our First Approach (‘Ea’ Representing the Results 
Obtained With Equation (16) and ‘Eb’ Representing the Results Obtained With Equation (18)) 

Fig. 6 depicts the results for the net electric field at P obtained using both equation (16) and 
equation (18). The graphs have been obtained using MATLAB R2014a. We note here that the 
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magnitudes of the net electric field for all the three values of cr 0 are considerably less (by one 

order of magnitude) than the corresponding values obtained without taking into account the 
effect of the backing plate. It is, thus, clear that the earthed backing disk has the effect of 
reducing the magnitude of the electric field produced at P. We also note that the difference 
between the results obtained using equations (16) and (18) is negligible for large distances of 
separation between the earthed sphere and the insulating disk. The difference starts becoming 
visible at about D = 80 mm for <r 0 = 1.5 x 10" 11 C/mm 2 , at about D = 40 mm for cr 0 = 0.5 x 

10' 11 C/mm 2 and at about D - 20 mm for cr 0 = 0.2 x 10 11 C/mm 2 . Also, the absolute 

difference is greater for greater values of the charge density. The results for the net electric 
field obtained from equation (18) are greater than those obtained from equation (16) by a 
maximum (occurring at D = 0.8 mm) of 25.4 % for all the three values of cr 0 . This difference 
goes down to 9 % at D = 20 mm and to less than 5 % at D = 40 mm for all the three values of 


Our final expression for the net electric field at P given by equation (34) (which is by 
far the most correct expression for this electric field) will produce a curve which lies in 
between the two curves (one curve obtained using equation (16) and the other curve obtained 
using equation (18)) shown for each value of cr 0 in fig. 6. The graphs that have been obtained 

for the net electric field at P using this expression for the three values of cr 0 that we have 

considered here (0.2 x 10’ 11 C/mm 2 , 0.5 x 10' 11 C/mm 2 and 1.5 x 10 11 C/mm 2 ) are shown 
below (in Fig. 7, Fig. 8 and Fig. 9 respectively). These graphs have been obtained using the 
Wolfram Development Platform (https://develop.open.wolframcloud.com/app/) . 
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Max. charge density = 0.2x10 A (-11) C/sq. mm. 



Fig. 7. Plot of the Net Electric Field at P Given By Equation (34) (Denoted Here By Enet), Plots of Its 
Components, Erl and Er2, Erl Being the Component Due to the Charges on the Charged Surface of the 
Insulating Disk and Their Images Produced Due to Their Reflection on the Earthed Metal Disk and the Images 
of All these Charges Produced in the Earthed Metal Sphere, and Er2 Being the Component Due to the 
Reflection of the Image Charges Produced in the Earthed Metal Sphere on the Earthed Metal Disk, and the 
Reflection of the New Image Charges Produced in the Earthed Metal Disk Due to the Abovementioned 
Reflection on the Earthed Metal Sphere, and the Reflection of the New Image Charges Produced in the Earthed 
Metal Sphere Due to the Aforesaid Reflection on the Earthed Metal Disk, and Such Reflections Going On on the 
Two Metal Surfaces Ad Infinitum, and a Plot of the Net Electric Field at P Given By Equation (16) Obtained 
Going By the Suggestion of Heidelberg (Denoted Here by Eh) for cr 0 = 0.2 x 10" 11 C/mm 2 
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Max. charge density = 0.5x10*(-11) C/sq. mm. 



Fig. 8. Plot of the Net Electric Field at P Given By Equation (34) (Denoted Here By Enet), Plots of Its 
Components, Erl and Er2, Erl Being the Component Due to the Charges on the Charged Surface of the 
Insulating Disk and Their Images Produced Due to Their Reflection on the Earthed Metal Disk and the Images 
of All these Charges Produced in the Earthed Metal Sphere, and Er2 Being the Component Due to the 
Reflection of the Image Charges Produced in the Earthed Metal Sphere on the Earthed Metal Disk, and the 
Reflection of the New Image Charges Produced in the Earthed Metal Disk Due to the Abovementioned 
Reflection on the Earthed Metal Sphere, and the Reflection of the New Image Charges Produced in the Earthed 
Metal Sphere Due to the Aforesaid Reflection on the Earthed Metal Disk, and Such Reflections Going On on the 
Two Metal Surfaces Ad Infinitum, and a Plot of the Net Electric Field at P Given By Equation (16) Obtained 
Going By the Suggestion of Heidelberg (Denoted Here by Eh) for er 0 = 0.5 x 10" 11 C/mm 2 
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Max. charge density = 1.5x10 A (-11) C/sq. mm. 



Fig. 9. Plot of the Net Electric Field at P Given By Equation (34) (Denoted Here By Enet), Plots of Its 
Components, Erl and Er2, Erl Being the Component Due to the Charges on the Charged Surface of the 
Insulating Disk and Their Images Produced Due to Their Reflection on the Earthed Metal Disk and the Images 
of All these Charges Produced in the Earthed Metal Sphere, and Er2 Being the Component Due to the 
Reflection of the Image Charges Produced in the Earthed Metal Sphere on the Earthed Metal Disk, and the 
Reflection of the New Image Charges Produced in the Earthed Metal Disk Due to the Abovementioned 
Reflection on the Earthed Metal Sphere, and the Reflection of the New Image Charges Produced in the Earthed 
Metal Sphere Due to the Aforesaid Reflection on the Earthed Metal Disk, and Such Reflections Going On on the 
Two Metal Surfaces Ad Infinitum, and a Plot of the Net Electric Field at P Given By Equation (16) Obtained 
Going By the Suggestion of Heidelberg (Denoted Here by Eh) for er 0 = 1.5 x 10" 11 C/mm 2 

Now, we shall find out the breakdown strength of an oxygen-enriched atmosphere 
considering our system shown in fig. 2 using Bower’s criterion [16]. The criterion is as 
follows. 


A E, 


•breakdown _ breakdown _ aXEj 

2v 4v 2 


breakdown 


(43), 


where 


v is the ionization potential of the surrounding gas (we take 12.5 V for its value for the case 
of oxygen), 

X is the mean free path of an electron (we take 3.7 x 10 7 in. for its value), 
and a is the radius of the earthed metal sphere (which is 0.8 mm. for our case). 

Putting these values in equation (43), we get 

3.38 X 10 7 

e E breakdown = 3.2 x 10~ 5 E breakdown - 4.7 x 10~ 13 E breakdown ... (44) 
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Plotting the expressions on the left-hand side and the right-hand side of this equation against 
Ebreakdown (using Wolfram Alpha Pro - http://www.wolframalpha.com) , we get the following 
graphs (shown in Fig. 10). 



Fig. 10. Graphically Solving Bower’s Equation (Equation (44)) Using Wolfram Alpha Pro 

( http://www.wolframalpha.com) 


We find in Fig. 10 that the graphs intersect at about E breakdown — 0.65 x 10 7 V/m. 
Thus, the solution of equation (44) is 


breakdown 


= 0.65 xlO 7 V/m. 


= 65 KV/cm. 

= 6.5 KV/mm. 

It may be observed from Fig. 9 that considering the value of the breakdown electric 
field of the oxygen gas in the regulator to be 6.5 kV/mm., an electrostatic discharge would 
ensue when the earthed metal sphere comes at a distance of (7.5 - 0.8) mm., i.e., 6.7 mm. 
from charged surface of the insulating disk if the maximum charge density ( a 0 ) on the disk is 

1.5 x 10" C/mm. . It thus appears that the occurrence of an electrostatic discharge inside an 
oxygen regulator is a possibility. 

Now, we can also calculate the minimum value of cr 0 which is required for a 
discharge to occur at point P using equation (42) and equation (41). For the dimensions we 
have used for our calculations, we have 

f 2 (S,a,b,R,s) = 8.675 [using equation (41)]. 
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Then, using s = 8.858 x 10 12 F/m = 8.858 x 10 15 F/mm and E breakdown = 6.5 KV/mm, we get 
<T) = [equation (42)] 

or,( j 0 =1.33x10 U C/mm 2 

This is the minimum value of cr 0 which is required for a discharge to occur at point P, the 

point of the earthed sphere nearest to the charged surface of the insulator and which lies on 
the perpendicular to the insulating disk at its centre. We may also call this charge density on 
the insulating surface the ‘critical charge density’. 

The minimum value of cr 0 which is required for a discharge to occur at point P when 
the centre of the earthed metal sphere is at a distance D from the charged surface of the 
insulating disk has been calculated for various values of D for the dimensions we have used 
here using equation (42 A). The graph thus obtained is shown in the fig. 11 given below. 



100 


200 


Din mm 


Fig. 11. The Minimum Value of cr 0 Which is Required for a Discharge to Occur at Point P When the Centre of 
the Earthed Metal Sphere Is at a Distance D from the Charged Surface of the Insulating Disk 

From the above results and discussions, we see that there is a possibility that an 
electrostatic discharge could occur inside an oxygen regulator. Now, considering the fact that 
the atmosphere inside an oxygen regulator is oxygen enriched for the most part, it also seems 
likely that such a discharge could also be responsible for fire hazards. 

The electric field inside the insulating disk has been calculated using equation (71) for 
the dimensions used here for the three values of cr 0 that we have considered here (0.2 x 10' 11 
C/mm 2 , 0.5 x 10 11 C/mm 2 and 1.5 x 10" 11 C/mm 2 ) and also for cr 0 = 1.33 x 10 ~ 11 C/mm 2 
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which has been found to be the critical charge density for various values of D and z. In doing 
these calculations, the relative permittivity of neoprene, which is 6.7 [17], has been used as 
the diaphragms of many oxygen regulators are made of neoprene [18]. Graphs for this 
electric field have been plotted using the Wolfram Development Platform 
(https://develop.open.wolframcloud.com/app/) . These graphs are shown below. 


Max. Charge Density oO = 0.2 x 10 A (-11) C/sq. mm., D = 0.8 mm. 

v 



- Erl 

- Er2 

— Enet = Erl + Er2 


Fig. 12. Plots of the Net Electric Field at a point z on the z-axis Inside the Insulating Disk (Denoted Here By 
Enet) Versus z and of Its Components, Erl and Er2, Erl Being the Component Due to the Charges on the 
Charged Surface of the Insulating Disk, i.e., the Surface of Charge B and Their Images Produced Due to Their 
Reflection on the Earthed Metal Disk, i.e., the Surface of Charge B' and Er2 Being the Component Due to the 
Image Charges Produced in the Earthed Metal Sphere Due to Reflection of the Charges on the Surfaces of 
Charge B and B' and the Reflection of the Image Charges Produced in the Earthed Metal Sphere on the Earthed 
Metal Disk, and the Reflection of the New Image Charges Produced in the Earthed Metal Disk Due to the 
Abovementioned Reflection on the Earthed Metal Sphere, and the Reflection of the New Image Charges 
Produced in the Earthed Metal Sphere Due to the Aforesaid Reflection on the Earthed Metal Disk, and Such 
Reflections Going On on the Two Metal Surfaces Ad Infinitum for cr 0 = 0.2 X 10 -11 C/mm 2 
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Max. Charge Density oO = 0.5 x 1 0 A (-11) C/sq. mm., D = 0.8 mm. 


— trl 

— Er2 

— Enet = Erl + Er2 



Fig. 13. Plots of the Net Electric Field at a point z on the z-axis Inside the Insulating Disk (Denoted Here By 
Enet) Versus z and of Its Components, Erl and Er2, Erl Being the Component Due to the Charges on the 
Charged Surface of the Insulating Disk, i.e., the Surface of Charge B and Their Images Produced Due to Their 
Reflection on the Earthed Metal Disk, i.e., the Surface of Charge B' and Er2 Being the Component Due to the 
Image Charges Produced in the Earthed Metal Sphere Due to Reflection of the Charges on the Surfaces of 
Charge B and B' and the Reflection of the Image Charges Produced in the Earthed Metal Sphere on the Earthed 
Metal Disk, and the Reflection of the New Image Charges Produced in the Earthed Metal Disk Due to the 
Abovementioned Reflection on the Earthed Metal Sphere, and the Reflection of the New Image Charges 
Produced in the Earthed Metal Sphere Due to the Aforesaid Reflection on the Earthed Metal Disk, and Such 
Reflections Going On on the Two Metal Surfaces Ad Infinitum for er 0 = 0.5 X 10 -11 C/?nm 2 
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Max. Charge Density 00 = 1.5 x 10*(-11) C/sq. mm., D = 0.8 mm. 



Fig. 14. Plots of the Net Electric Field at a point z on the z-axis Inside the Insulating Disk (Denoted Here By 
Enet) Versus z and of Its Components, Erl and Er2, Erl Being the Component Due to the Charges on the 
Charged Surface of the Insulating Disk, i.e., the Surface of Charge B and Their Images Produced Due to Their 
Reflection on the Earthed Metal Disk, i.e., the Surface of Charge B' and Er2 Being the Component Due to the 
Image Charges Produced in the Earthed Metal Sphere Due to Reflection of the Charges on the Surfaces of 
Charge B and B' and the Reflection of the Image Charges Produced in the Earthed Metal Sphere on the Earthed 
Metal Disk, and the Reflection of the New Image Charges Produced in the Earthed Metal Disk Due to the 
Abovementioned Reflection on the Earthed Metal Sphere, and the Reflection of the New Image Charges 
Produced in the Earthed Metal Sphere Due to the Aforesaid Reflection on the Earthed Metal Disk, and Such 
Reflections Going On on the Two Metal Surfaces Ad Infinitum for <r 0 = 1.5 X 10 _11 C/mm 2 
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Max. Charge Density 00 = 1.33 x 10*(-11) C/sq. mm., D = 0.8 



Fig. 15. Plots of the Net Electric Field at a point z on the z-axis Inside the Insulating Disk (Denoted Here By 
Enet) Versus z and of Its Components, Erl and Er2, Erl Being the Component Due to the Charges on the 
Charged Surface of the Insulating Disk, i.e., the Surface of Charge B and Their Images Produced Due to Their 
Reflection on the Earthed Metal Disk, i.e., the Surface of Charge B', and Er2 Being the Component Due to the 
Image Charges Produced in the Earthed Metal Sphere Due to Reflection of the Charges on the Surfaces of 
Charge B and B' and the Reflection of the Image Charges Produced in the Earthed Metal Sphere on the Earthed 
Metal Disk, and the Reflection of the New Image Charges Produced in the Earthed Metal Disk Due to the 
Abovementioned Reflection on the Earthed Metal Sphere, and the Reflection of the New Image Charges 
Produced in the Earthed Metal Sphere Due to the Aforesaid Reflection on the Earthed Metal Disk, and Such 
Reflections Going On on the Two Metal Surfaces Ad Infinitum for er 0 = 1.33 X 10 -11 C/mm 2 

The breakdown electric field strength of neoprene lies between 15.7 MV/m. to 26.7 MV/m. 
[19], i.e., between 15.7 KV/mm. to 26.7 KV/mm.. We can see from the graphs above that the 
electric field inside the neoprene diaphragm does not reach such values for values of a 0 upto 
1.5 x 10" C/mm". The electric field inside the neoprene diaphragm is expected to reach the 
breakdown electric field strength of neoprene at much higher values of er 0 . 


7. Conclusion 

Thus, in trying to investigate the electric field which is likely to be produced at the tip of the 
valve stem inside a typical oxygen regulator due to the charge that might be produced on the 
diaphragm, we have first modelled the system as a very small earthed metal sphere in front of 
a charged insulating disk partially backed by an earthed metal disk. We have qualitatively 
argued that the density of charge that could be produced on the diaphragm might suitably be 
represented by the Gaussian function. We have then derived an analytical expression 
(equation (34)) for the net electric field produced at the point on the earthed metal sphere 
nearest to the charged surface of the insulating disk. Next, we have calculated the magnitude 
of the electric field produced at this tip of the earthed sphere in our model using 
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measurements obtained from a sample oxygen regulator for the dimensions of the various 
components of our model. We have also found out an expression for the critical charge 
density on the surface of the insulating disk (equation (42)), and have also found out a value 
for this critical charge density for the sample oxygen regulator we have considered, which is 
1.33x10 11 C/mm 2 . We have also found out an expression for the electric field inside the 
insulating disk in our system (equation (71)), and have calculated values of it for the sample 
oxygen regulator we have considered here. It appears from our analyses that there is a 
possibility that an electrostatic discharge could occur inside an oxygen regulator, and due to 
the presence of an enriched-oxygen atmosphere, such a discharge might also lead to a fire 
incident. 
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